function f = bisect_wage(y,W0,U0)

x = exp(y) ;

f = ((W0(x)-U0).^2) ;